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Motivated by the recent experiment of Ref. [1] , we develop a mean- field theory of the interaction- 
induced antiferromagnetic (AF) state in bilayer graphene at charge neutrality point at arbitrary 
perpendicular magnetic field B. We demonstrate that the AF state can persist at all B. At higher 
B, the state continuously crosses over to the AF phase of the f = quantum Hall ferromagnet, 
recently argued to be realized in the insulating i/ = state. The mean-field quasiparticle gap is finite 
at _B = and grows with increasing B, becoming quasi-linear in the quantum Hall regime, in accord 
with the reported behavior of the transport gap. By adjusting the two free parameters of the model, 
we obtain a simultaneous quantitative agreement between the experimental and theoretical values 
of the key parameters of the gap dependence - its zero-field value and slope at higher fields. Our 
findings suggest that the insulating state observed in bilayer graphene in Ref. is antiferromagnetic 
(canted, once the Zeeman effect is taken into account) at all magnetic fields. 

PACS numbers: 



Introduction. Bilayer graphene (BLG) presents an ex- 
citing arena for the observation of the correlated electron 
physics [Tk23j. Nearly quadratic dispersion of the elec- 
tron spectrum about the charge neutrality point makes 
the system susceptible to even weak interactions and at 
zero magnetic field, allowing for instabilities towards var- 
ious broken-symmetry phases. A variety of correlated 
states at zero doping, characterized by different order- 
ing of the valley, layer, and spin degrees of freedom have 
been predicted or considered ITS] . At finite perpen- 

dicular magnetic field B, quenching of the kinetic energy 
facilitates the correlation effects. In the quantum Hall 
(QH) regime, the zero-density state transforms into the 
ly ^0 quantum Hall ferromagnet (QHFM) [THU HMS] , 
which also supports a number of interesting phases. 

Recent transport experiments [H [SHZ] on high-quality 
suspended BLG samples provided compelling evidence 
for the interaction-induced ground states both at _B = 




FIG. 1: The antiferromagnetic (AF) state in bilayer graphene 
(BLG) at arbitrary orbital magnetic field. If the Zeeman ef- 
fect is neglected, as done in this paper for simplicity, the A 
and B sublattices, located in different layers, have arbitrary 
antiparallel spin polarizations, as shown. Once the Zeeman 
effect is included, the AF state transforms into the canted AF 
state p3], not shown here. 



and in the QH regime. Several qualitatively different 
behaviors were reported. In Refs. [HIISI, the zero-density 
state was insulating in the QH regime (reached already at 
B > IT), showed metallic value of the two-terminal con- 
ductance G > /h a,t B = and a nonmonotonic behav- 
ior of G at intermediate B < IT. In Ref. [7], in cleaner 
samples (labeled B2 therein), at _B = 0, the differential 
conductance displayed signatures of the insulating gap 
with the minimal zero-bias conductance G « 0.2e^/h; 
remarkably, at the same time, no fully developed insu- 
lating state was observed at higher B. Finally, Ref. [1] 
reported a pronounced insulating state at all magnetic 
fields. The transport gap was £'™p w 20K at B = 
and grew with increasing B, becoming linear in B in the 
QH regime, with the slope dE^^^/dB « 5.5meV/T, as 
the state continuously crossed over to the v = QHFM 
state. 

While all scenarios are equally interesting, in this pa- 
per we concentrate on the theoretical description of the 
latter [T] - the insulating state at all magnetic fields. We 
develop a mean-field (MF) theory of the insulating an- 
tiferromagnetic (AF) state in BLG, Fig. [l] at arbitrary 
perpendicular magnetic field. We demonstrate that, the 
AF phase can persist at all B, continuously interpolating 
between the earlier studied B = O HSl IIIl |16| and QH 
(AF phase of the i/ = QHFM) [13] limits. Most impor- 
tantly, the obtained mean-field spectrum reproduces well 
the crucial experimental feature of Ref. pTj - the depen- 
dence of the transport gap on the magnetic field. We ob- 
tain a simultaneous quantitative agreement between the 
experimental and theoretical values of the key parame- 
ters of the gap dependence - its zero-field value and slope 
at higher fields - by adjusting the two free parameters of 
the model. Our findings further substantiate the conclu- 
sions about the AF phase in the QH regime [23] and at 
S = [T], suggesting the AF phase as the most likely 
candidate for the insulating state observed in Ref. [T]. 
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AF phase in BLG. A large number of correlated phases 
in BLG at -B = predicted or considered in theoretical 
literature [SHIBj [18] can be classified according to the 
properties of their charge excitations as (i) bulk gapless 
(e.g., nematic [SJ [TH [131 HH] ) ; (ii) topologically nontrivial 
bulk gapped phases with gapless edge excitations (e.g., 
quantum anomalous Hall (QAH) [T^HTOl [2T] , quantum 
spin Hall (QSH) [15l[16l[18]); and (iii) fully gapped (bulk 
and edge) (e.g., ferroelectric [H [11 [H] , AF |l [H [TH 
[Tni[TH]). The AF phase was argued in Ref. [T] to be the 
most likely candidate for the insulating state at B = 0, 
while the phases of (i) and (ii) types can be ruled out with 
certain confidence, since they should exhibit metallic two- 
terminal conductance G > e^/h. The phases (i) or (ii) 
are more suitable candidates for the metallic low-field 
behavior observed in Refs. [51 [B] . 

In the QH regime, the zero-density state transforms 
into the 1/ = QHFM [HHSS]. The generic phase di- 
agram of the 1^ = QHFM in BLG was obtained in 
Ref. [23] and consists of four phases: spin-polarized, an- 
tiferromagnetic (canted, once the Zeeman effect is taken 
into account), interlayer-coherent (at zero perpendicu- 
lar electric field), and fully layer-polarized. Also, it was 
argued in Ref. [23] that the experimentally observed in- 
sulating ly = QH state in BLG is the AF phase of the 
V — QHFM. This conclusion was reached by compar- 
ing the obtained phase diagram with the experimental 
data of Ref. [5] and was based on the argument that 
AF is the only phase consistent with the observation of 
the insulator-insulator phase transitions in the perpen- 
dicular electric field. The same transitions are observed 
in Ref. [Ij and thus the same conclusion about the AF 
phase in the QH regime can be made. 

Crucially, combined with the above conclusions, the 
fact that the insulating state of Ref. |lj shows a contin- 
uous crossover between the zero-field and QH regimes 
strongly suggests that the AF phase persists at all mag- 
netic fields. Here we theoretically demonstrate that this 
is indeed a feasible scenario. 

Model. Our starting point is the Hamiltonian for in- 
teracting electrons in the perpendicular magnetic field, 

H = Hq + Hi, 

^0 = 1 dhij^ho^, ho = ^ {%+Pl+%Jl) ,(1) 

= ^ / E « ^ ■■[^^To.pi^r-- ■ (2) 

^ / cx-P Til 

We describe electron dynamics in the framework of the 
two-band model [24] of BLG, valid at energies e ^ t_L 
below the interlayer hopping amplitude ij^ w 0.3eV. At 
such energies, the wave-functions are predominantly lo- 
calized on A and B sublattices, located in different lay- 
ers. Fig. [l] The relevant degrees of freedom are joined 
into the eight-component field operator ^ = {tp^^ip^Y, 



product KK' (E) AB ig) s of the valley, sublattice, and 
spin spaces, respectively. We use the same basis as in 
Refs. [121 HH] ■ Note that in this basis, the actual A and 
B sublattices are interchanged in the K' valley; there- 
fore, to avoid confusion, we denote this sublattice space 
as AB. In Eq. ([2|, : . . . : denotes normal ordering of op- 
erators and the summation goes over a, P € {0, x, y, z}. 
In the kinetic energy term ([!]), m is the effective mass, 

'P± =Px± iPy, Pa =Pa~ %^a. Pa = -I^q for a = X,y, 

and rot A = (0, 0, B). In Eqs. and ^ and below, for 
a,/3,7 € {0,x,y,z}, 

with the unity (tq = 1) and Pauli {tx, Ty, Tz) matrices 
in the corresponding subspaces. To keep the analysis 
simpler, we leave the orbital magnetic field as the only 
single-particle effect and neglect the effects of warping 
and strain [5S]: the quite large extracted value -Egap ~ 
20K of the transport gap suggests that the correlation 
effects dominate over these effects under the experimental 
conditions of Ref. [1] . We also neglect the Zeeman effect 
for the same reason: for perpendicular field orientation, 
the actual canted AF phase should differ little from 
the AF phase. 

Equation ^ is the most general form of the two- 
particle interactions, asymmetric in the KK' ®AB space, 
allowed by the symmetry of the BLG lattice [ni[T51[^ . 
The couplings satisfy the relations gj^j^ = g^x = 9xy = 
9yx — 9yyi 9-l-z = 9xz — 9yz} 9z± = 9zx — 9zy, 
9±a = 9x0 = gyo, 9q± = 9ox = goy, yielding the total of 
9 independent couplings [l2l [131 [HI ES] . The symmetric 
coupling goQ characterizes the strength of the screened 
Coulomb interactions [27]. For static screening in the 
large- approximation, goo = l/(2iVln4), with the ac- 
tual N = 4: due to two spin projections and two valleys. 
The asymmetric interactions can arise from both the ac- 
tual Coulomb interactions and electron-phonon interac- 
tions. 

There is no accurate knowledge of the coupling con- 
stants 9ai3i yet their "bare" values at the bandwidth 
^ t± of the two-band model determine the favored 
broken-symmetry ground state in BLG at zero doping. 
At B = 0, a systematic weak-coupling analysis of the 
many-body instabilities is carried out within the RG ap- 
proach [HI [121 [131 HH] . In the QH regime, the interaction- 
induced v = state is studied within the framework 
of QHFMism [T9H23] . Among the variety of predicted 
phases, the AF phase was demonstrated to occur at both 
B ~ and in the QH regime, under realistic assump- 
tions about the values of the coupling constants gafi . We 
will now assume that the AF phase is the favored ground 
state both at B = and in the QH regime and demon- 
strate that the AF phase then persists at all intermediate 
B and that the two limits are adiabatically connected. 

Mean-field analysis. We study the problem within the 
MF approach. At _B = 0, the order parameter (OP) 
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FIG. 2; The components Ao 



'es and of the mean- 



field potential S of the antiferromagnetic phase ( AF) as func- 
tions of es/Aj; Az is obtained by numerically solving the 
self-consistency equation, either Eq. (JSl or (|10||. The value 



g-~- = 0.4 was used. See caption to Fig. [3]for details 
Q = {-.ipilj^ :) of the AF phase has the form 



AB 



(3) 



At finite magnetic field B > 0, due to the emergence of 
the n = 0, 1 Landau levels (LLs) [24] and the peculiar 
property of their wave-functions to reside on only one 
sublattice in each valley, the OP necessarily acquires a 
component t^^ ~ ~ Tz)^^ in the AB space; hence 
one needs to include the TzOz component in the full OP 
of the AF phase. The OP that describes the AF state at 
arbitrary magnetic field therefore has the form 



Q = rf ^ <E) {Qol + QzTz) 



AB 



(4) 



Performing decoupling of interactions in Eq. ([2| , Hi 
Hi^mi = J d^rtp'^ Alp, we obtain the MF potential 



.KK' 



(Aoi 



\AB 



(5) 



where A„ = - ^g''°"'Qa [a = 0, z) and g"^"^^ = goo+gzz + 
45_L_L - 2g_Lz - 2gz± - 2go± + goz - 2g±o + gzo and g''^^ = 
goo + gzz - 'ig±± - 2g±z + 2gz± + '2go± + goz - 2g±Q + g^o- 
Solving the eigenvalue problem for the Hamiltonian 
ho + A, we obtain the mean-field LL spectrum 



EnXa = (Ao - A^)sAScr, n = 0,l, 

E,,±xa = AqSas. ± V4 + A2, n > 2. 



(6) 
(7) 



Here, e„ = e_B-\/ n{n — 1) is the LL spectrum of the non- 
interacting BLG and ejs = l/{ml%) « 1.3^^B[T]K is 
the cyclotron energy , in which Is — \/ rnc/ {eB) is the 
magnetic length and rrie is the electron mass. Each state 
is characterized by the valley X — K,K' and spin a =t, i 
indices and = ±1 and Sg- = ±1, respectively. 

Calculating the OP ^ in the eigenstate basis, we ob- 
tain the self-consistency equations at zero temperature 



A, 



A, 



g 



1 



Ao = 



-g es- 



(8) 

In the right-hand side of the equation for A^, the unity 
represents the contribution from n = 0, 1 LLs, while the 
sum from n > 2 LLs. We impose an ultraviolet energy 




FIG. 3: Mean-field Landau level (LL) spectrum (|6|,([7| of the 
antiferromagnetic (AF) phase as a function of the magnetic 
field B, obtained using the numerical solution for Ao and A^ 
shown in Fig. [2] the dependence on B is expressed in terms 
of the ratio e_g/A° of the cyclotron energy es = eB/{mc) 
and zero-field gap A*^. Only the levels with n < 4 are 
shown (the numbers indicate n), E„±Kt = En±K'i (red), 
E„±Ki = En±K'^ (blue). The dependence of the gap -Egap = 
2 min(_E2+ifti E'Oift) between the positive and negative en- 
ergy levels on es closely reproduces that of the transport gap 
£'gap of the insulating state observed in Ref. [T], compare with 
Fig. 3 therein. The used value g"^^ = 0.4 provides quantita- 
tive agreement with the experimental value dE^^^ / des = 1.3 
(converted from d£™P/dB = 5.5meV/T at m = 0.028me [6]) 
of Ref. [I ] for the slope diSgap/des of the gap at higher fields, 
see Eqs. <\13\, (14 1 and text. 



cutoff eo on the spectrum and cut the otherwise log- 
divergent sum by the large integer part no — [eo/es] ^ 1. 

Eqs. ([4|-([8| are the key result of our work. The solu- 
tion of Eq. Isl) for Az and Aq determines the evolution 
of the order parameter (|4| of the AF state with the mag- 
netic field and the MF quasiparticle spectrum (|6|, ([7|. 
Below we discuss the key properties. 

At i? = 0, Ao = and the equation for A^ reduces to 



g 



de 



A, 



v/^^+Al 



(9) 



Its solution A° — 2eo exp(— 1/g^^^) determines the gap 
in the spectrum at B = 0. The OP is given by Eq. ^ 
with Q° = -A°m/{AT:g^^^). 

One can eliminate the cutoff no from Eq. ([s]) for Az 
using the standard procedure known from the ECS the- 
ory |28| . Namely, one may represent the integral 



de 



y^^+Al 



"0 

n=2 



In 



eg 
A. 



o(l) 



by an arbitrary series with asymptotic /o + '^^"=2 /« ~ 
ln(2no) -t- o(l) at no = [eo/es] — > oo; the specific form 
of the series is a matter of convenience. Adding and 
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subtracting these two forms from the sum in Eq. ([8]) , one 
arrives at an equivalent equation for , 




In 



1 \ . es 

A?- 

(10) 

This form shows exphcitly that the functional depen- 
dence of Az on es is, in fact, fully determined by one 
parameter, its value A° at B = 0. The component 
Aq = — g^^^EB, in its turn, is linear in B and its slope 
is controlled by the coupling constant (?^°^. The present 
theory is therefore described by two parameters, the zero- 
field gap A° and the coupling constant g^^^ 



In the QH regime es ^ A^, it follows from Eq. (10) 



A, 



r'{eB)eB, .r"'(eB) = l/[ln(es/A°) + Co], (11) 



Co = H'^e^'^)+J2n=2[^/n-l/Vr>'in-l)] » 0.674. We 
have introduced the notation ^^^^(es), since according to 
Eq. (|8]), Eq. (Ill can be interpreted as a QH-type depen- 
dence, akin Ag — — g^'^^es, with a renormalized coupling 
constant 5^^^ ^ ^^^^(£5) = g^^V[l-g^^^ ln(1.02eo/eB)]. 

Thus, at higher fields, A^ is quasi-linear in cb with a 
logarithmically varying slope. The OP equals 



1 



5g 



zzz g'^^'^ 



9^ 



(12) 

The first contribution arises from 



n = 0, 1 LLs, while the second one is the AF OP induced 
in the n > 2 LLs by LL mixing. This OP describes the 
AF phase of the z/ QHFM considered in Ref. [5^. 
At intermediate fields, either Eq. (Isl) or (10) for A^ 



can be solved numerically. The components Aq and A^ 
as functions of the magnetic field B expressed in terms 
of es/A^ are plotted in Fig. [2] As anticipated, we find 
that upon applying the magnetic field the system pre- 
serves the AF order and the AF state at i? = [Eq. ^] 
continuously crosses over to the AF phase of the v — Q 
QHFM [Eq. 

The resulting mean- field LL spectrum (|6]),([7| of the 
AF state is plotted in Fig. |3] The MF quasiparticle gap 
Sgap = 2T[mi{E2+K\,EDK\) is given by twice the energy 
of the lowest positive state. The gap E'gap = 2A" is 
finite at B = 0. As seen from Figs. [2] and |3j apart from 
a region es < A^ around _B = 0, in the major range of 
fields e_B ^ A°, the gap is determined by n = 0, 1 LLs, 

i?gap/2 = iJoA'T = A. + IAol « [r"^(eB)+.g""1eB, (13) 



and the formula ( 11 ) is accurate. Thus, in the QH regime. 



reached already at es > A^, the gap (13) has a quasi- 



linear dependence on associated with the QHFM 
physics; the quasi-slope equals 



d£^gap ^2 '^^^^^ 



de. 



de. 



9 



dA, 
des 



(14) 



Comparison with the experiment of Ref. \T^. The ob- 
tained spectrum. Fig. [3j reproduces well the dependence 
of the transport gap E'^{B) on the magnetic field B 
observed in Ref. [T], compare with Fig. 3 therein. A dis- 
tinct feature of E'^{B) is that the linear dependence at 
higher S, if extrapolated to S = 0, crosses the vertical 
axis at a value only slightly below the actual zero-field 
gap. The obtained spectrum in Fig. [3] exhibits the same 
property. 

A quantitative agreement between the key parame- 
ters of the experimental E'^{B) [T] and theoretical 
-£'gap(-B) (Fig. [3]) gap dependencies, the zero- field value 
E^^IIb = 0) « 2GK and the slope AE°^l/deB = 1.3 at 
higher fields |29j, is achieved by adjusting the two free 
parameters of the model, A° and g^^^ . First, the zero- 
field gap is fit by setting A^ = £;|^p(B = 0)/2 w lOK. 
Second, taking the typical slope dA^/des ~ 0.25 of the 
A^(eB)-dependence in the experimentally relevant range 



IT < B < 4T of fields, we obtain from Eq. (|l4j) that the 
experimental slope is fit at « 0.4. The Az{tB) and 
Ao(eB) dependencies in Fig.[2]and the spectrum in Fig.js] 
are presented for this value of 3^°^. 

In this regard, we mention that after a preprint |30| 
of the present work became available, a similar MF 
analysis of the AF phase at finite B was performed 
in Refs. [3T1 However, in contrast to the present 

study, no quantitative agreement between the theoret- 
ical value of the slope at higher fields and the exper- 
imental value dE^^l/AiB = 1-3 of Ref. [I] could be 
achieved there. This discrepancy originates from disre- 
garding the Aq component of the MF potential ([5| of the 
AF phase in the analysis of Refs. [SI] |32] and consider- 
ing only the A^ component. According to the approach 
of Refs. [SD ISa the gap E^Y/^ = 2A^ of the AF state 
is determined solely by A^; then its maximum possible 
slope max(d£'^p^/deB) — 0.5 is indeed smaller than the 
experimental value by a factor of 2.6. We emphasize 
that the component Aq is automatically generated by 
the wave-functions of the n = 0, 1 LLs and without it 
the self-consistency gap equations ([s]) simply cannot be 
satisfied. Once Aq = —g^^^eB has been taken into ac- 
count, both components Aq and A^ contribute to the 



gap ( 13 ) and its slope ( 14 1 can always be fit by adjusting 
the coupling as we demonstrated above. 

Conclusion. We developed a mean-field theory of the 
interaction-induced antiferromagnetic state in BLG at 
charge neutrality point at arbitrary perpendicular mag- 
netic field. The theory reproduces well the key features 
of the recent experiment [T] on suspended BLG samples: 
persistence of the insulating state at all magnetic fields 
and the dependence of its transport gap on the magnetic 
field. At higher magnetic fields, the state crosses over to 
the antiferromagnetic phase of the 1^ = QHFM, argued 
in Ref. |23| to be realized in the insulating v = quan- 
tum Hall state. The presented analysis suggests that the 
insulating state observed in Ref. T is antiferromagnetic 
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(canted, once the Zeeman effect is taken into account) at 
all magnetic fields. 
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